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This paper presents a state estimation approach for an uncertain linear equation with a non-invertible operator in Hilbert space. 
The approach addresses linear equations with uncertain deterministic input and noise in the measurements, which belong to 
a given convex closed bounded set. A new notion of a minimax observable subspace is introduced. By means of the presented 
approach, new equations describing the dynamics of a minimax recursive estimator for discrete-time non-causal differential- 
algebraic equations (DAEs) are presented. For the case of regular DAEs it is proved that the estimator's equation coincides 
with the equation describing the seminal Kalman filter. The properties of the estimator are illustrated by a numerical example. 
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1 Introduction 

The importance of models described by DAEs (or de- 
scriptor systems) in economics, demography, mechanics 
and engineering is well known [13]. Here, motivated 
by further applications to linear DAEs, we present 
a state estimation approach for linear deterministic 
models described by an abstract linear equation in a 
Hilbert space. Our approach is based on ideas under- 
lying H^jH^ filtering [4,3] and set-membership state 
estimation [7,14,16,12,15]. 

-/^-estimators like Kalman or Wiener filters [4,1] give 
estimations of the system state with minimum error 
variance. The ^-estimation problem for linear time- 
variant DAEs was studied in [17] without restricting the 
DAE's matrices. The resulting algorithm requires the 
calculation of the so-called "3-block matrix pseudoin- 
verse" . In [10] the authors introduced explicit formulas 
for the 3-block matrix pseudoinverse and derived a 
recurrence filter, assuming a special structure for the 
DAE's matrices. A brief overview of steady-state H 2 - 
cstimators is presented in [8]. 

Hoc estimators minimize a norm of the operator map- 
ping unknown disturbances with finite energy to filtered 
errors [3]. We stress that the Hoo estimator coincides 
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with a certain Krein space H 2 filter [19]. The H^ fil- 
tering technique was applied to linear time-invariant 
DAEs with regular matrix pencils in [22] . 
A basic notion in the theory of set-membership state 
estimation is that of an a posteriori set or informational 
set. This notion has roots in control theory [4]. By def- 
inition, it is the set of all possible state vectors p, that 
are consistent with a measured output y, provided that 
an uncertain input / and measurement error rj belong 
to some bounded set Sf . We will be interested in the case 
when the state ip £ H. obeys an abstract linear equation 
Lp = /, provided y = Up + ry, (/, rj) G <S ', where Sf is 
a bounded closed convex subset of an abstract Hilbert 
space. The problem is to find an estimation <p of p with 
minimal worst-case error. This problem was previously 
considered in [16,14]. Due to [16] a vector p is called 
a linear minimax a-posteriori estimation (or a central 
algorithm due to [14]) iff W e rl 



(£,0) 



(sup (£,<p) 



■w£U,<p))/2, 



(1) 



provided that an a posteriori set &(y) := {p : (Lp,y — 
Hip) <G Sf } is a bounded convex subset of the Hilbert 
space rl. Note that if there exists po so that Lp = 0, 
Htpo — 0, then sup^( y \ (£, p) = +oo for some L Thus, the 
above approach does not work if L is non-injective. In 
this paper we generalize the approach of [14,16] to linear 
equations with non-injective L. Futher generalization is 
presented in [25]. 
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The main contribution of this paper is a new notion of 
a minimax observable subspace C for the pair (L, H ) 
(Definition 1). It is useful when one needs to evaluate 
a priori how far the estimation tp is from a "real" state 
ip in the worst case, provided tp is constructed from the 
measurements y. Due to Proposition 1, the worst-case 
estimation error is finite iff C — H; otherwise tp may be 
too far from a "real" state tp for some directions £ £ H, 
even for bounded / and rj. In fact, given y, we can pro- 
vide an estimation with finite worst-case error for the 
projection of tp onto C only. Thus £ is an analog of the 
observable subspace [4, p. 240] for the pair (L, H) in the 
context of set-membership state estimation. 
The introduced notion allows the generalization of ideas 
from [14,16,21] to non-injective linear mappings, in par- 
ticular for the case Lp(t) — {{Ftp) t — C(t)p(t),Fp(t )) 
with F £ K mxn which arise in the state estimation for 
linear continuous non-causal : DAEs [24] . As a conse- 
quence, one can apply the minimax framework, orig- 
inally developed [7,5] for DAEs (F = E in the linear 
case) with bounded uncertainties, to DAEs [25] with 
unbounded inputs (see example in Section 3). 
In order to stress connections with H x approach, we 
note that the minimax framework [7,5] incorporates the 
set-membership state estimation and H^ filtering for 
ODEs by application of dynamic programming [2] to the 
informational state X(t): for linear ODE the worst-case 
estimation is set to be the Tchebysheff center of X(t). 
Although we derive the estimation from the minimiza- 
tion of the worst-case error, as it is stated in Definition 
1, our approach (for the ellipsoidal bounding set §f and 
causal 2 DAEs) results in the same estimation and error 
as in [7]. Thus, the f-minimax estimation gives a proper 
generalization of the recurrence algorithm from [7] to 
the case when X(t) may be unbounded. 
We illustrate the benefits of the new notion by intro- 
ducing a minimax recursive estimator for discrete-time 
non-causal DAEs: it works for non-causal DAEs un- 
like [10,8,22,23] and for the regular case it coincides 
with one proposed in [10] (Corollary 1). In addition, the 
minimax observability subspace allows one to identify 
the observable (in the minimax sense) part of the state 
with respect to given measurements. Computing the 
index of non-causality, one can a priori check how good 
connections between observations and state are: models 
with zero index are fully observable while models with 
non-zero index have an unobservable part in the state. 
This paper is organized as follows. In Section 2 we give 
definitions (Definition 1) of the minimax estimation, er- 
ror and observable subspace for abstract linear equations 
and we construct the estimation for a convex bounded 
Sf, in particular for an ellipsoidal Sf (Proposition 1). 
In Section 3 we introduce the minimax observable sub- 

DAE is said to be non-causal if the corresponding initial- 
value problem has more than one solution. 
2 Note that the dynamic programming was previously ap- 
plied to causal DAEs in [6] in order to construct a regulator 
in LQ-control problems with DAE constraints. 



space and index of non-causality for DAEs in discrete 
time (Definition 2) and we derive the minimax estimator 
(Theorem 1). Also we discuss connections to FLijFL^ 
framework (Corollary 1) and present an example. 
Notation. Linear mappings: (■,■) denotes the inner 
product; Jf("Hi,^ 2 ) denotes the space of all bounded 
linear mappings from Hi to H2, -£?(%) := Jz?{H 7 H); 
t-H is the identity mapping in H; 2>{L), R(L), N(L) 
denote, respectively, the domain, range, and null-space 
of a linear mapping L : @(L) n- R(L); L* : H —>■ H 
is the adjoint of L; F' denotes the transpose of F; F + 
is the pseudoinverse of F; E is the identity matrix; 
diag(^4i . . . A n ) denotes a diagonal matrix with A t , 
i = l,n on its diagonal; {a^s}™ := (x\, . . . , x n ) is an 
element of Hi x • • • x H n , 0„ m £ R mxn denotes the 
m x n-zero matrix. 

Functionals: I\(x) := (QiLx,Lx) + (Q2Hx,Hx), 
I(x) := (QiLx, Lx) + (Q2(y — Hx),y — Hx) with pos- 
itive definite self-adjoint Qi £ -^{F) and Q2 £ -£?(30i 
c(G, x) := sup{(:r, y),y £ G} is a support function of G; 
>Y±:=l(c(9(3,),£)±cMy),-e));\\x\\ 2 s = (Sx,x). 
Sets:^ l3 {0) := {x : h(x) < /3},Sf(0) := {x : h{x) < 1}; 
dom/ := {x : f(x) < 00} is an effective domain of /; 
Arguing/ := {x : f(x) = min x /} is the set of global 
minima of /; G is the closure of a set G. 



2 Linear minimax estimation problem in a 
Hilbert space 

Let vector y be observed in the form of 



y = Hp + r] 



(2) 



where p obeys the equation 



Ltp = f (3) 

We assume that L : 3>(L) H> T is a closed linear op- 
erator [4, p. 63], @{L) C H is a linear set, @(L) = H, 
H <G J2?(H, y). Also we assume that (/, rf) is an unknown 
clement of a convex bounded closed set Sf C F x y, 
H, F, y are Hilbert spaces. 



Definition 1 Let C 



E H : p(£) < 00} with 



p(£)~ igf p(t,<P), P&<P):= sup \(£,tp-i>)\ 

f<^^{y) tpe^(y) 

The set C is called a minimax observable subspace for 
the pair (L, H) . A vector tp £ &(y) is called a minimax a 
posteriori estimation in the direction^ £ (£- estimation) 
if p(£, tp) = p(£). The number p(£) is called a minimax a 
posteriori error in the direction £ (£-error). 



3 The subspace I := {a£ : a £ R}, assigned with £ £ H, 
defines some direction in H. To estimate tp in the direction 
£ means to estimate the projection (£, tp)£ of tp onto L 



Our aim here is, given y, to construct the ^-estimation 
p of the state <p, terror p{£) and minimax observable 
subspace C, provided t £ C. Note that p{£) = +00 if 
£ £ C so that any ip £ H. is a ^-estimation by Definition 1. 



Proposition 1 Let Sf he a convex closed bounded subset 
o/JxJ. Then t £ C iff I, -I £ domc(S%), •)• If £ £ C 
then the (.-estimation (p along with t - error p{€) obey 



(£,0) = 1 _,p(£) = 1+ 



(4) 



Define T : ®{T) ^ H by the rule T(z, u) := L* z + H*u 
with &(T) := 9{L*) x y and let 

^ = {(f,v)-(Qif,f) + (Q2V,v)<n 



If R(T) = R(T) then x £ Argmin^ / is the (.-estimation, 
C = domc(Sf(0), •) = R(T) and 

p(£) = (1-I(x))*c(&(0),£). (5) 

The worst-case estimation error for any direction is 



(6) 



sup \\x — x\\ = inf sup \\(p — x\\ = 
(1-J(x))2 sup c(&(0),£) = sup p{£) < +00 



If C = % then (6) is finite. 



PROOF. Let £(S%)) = {(£,ip),ip£ S%)}. Since ^(y) 
is convex (due to convexity of 2>(L) and Sf) and x 1— > 

(£, x) is continuous, it follows that £{fS(y)) is connected. 
Noting that xai<g< v \(l, ip) = —c(&(y), —£), we see 



£(<%)) = [-c(^(y),-£),c^(y) 1 £)]cR i 

Thus p(£, ip) = +00 if £, -£ £ dom c(&(y), •)• Otherwise 
l(@(y)) is bounded, implying £ £ C. Hence, £ £ L iff 

Let £ £ £, p £ @(y)- Since £(^(y)) is connected, there 
exists y* € ^(y) so that (£, y>„) =7- is the central point 
of £(@(y)). The worst-case distance p(£,(p) is equal to 
the sum of the distance \(£,p> — <p*)\ between (£, tp) and 
the central point 7_ and the distance 7+ between one of 
the boundary points of l(fS{y)) and 7_. Therefore, 7_ 
has the minimal worst-case distance 7+. Hence, p>„ = p 
due to Definition 1, which implies (4). 
We proceed with the ellipsoidal &. Let Q\ = Ijr, Q 2 = 
t y for a simplici ty. Due to [11, Sec 5.§3] R(T) = R(T) 
implies R(T*) = R(T*). Thus [4, p.l4,Cor.l.4.3], there 
exists x £ S'(L) so that T*x is the projection of (0,y) 



onto R(T*) = {(Lx,Hx),x £ 3>(L)}, implying x £ 
Arguing/, and 

(y - Hx, Hx) = (Lx, Lx) , Vx £ 3>{L) 

Noting this, one easily derives 4 I(x — x) = 1^ (x) + I(x) 
for all x £ @{L). Having it in mind and noting that 
@(y) = {p : I{p) ^ 1} and x £ @(y), one derives 

£ + gf^(0)=Sf(y) (7) 

where /3 := 1 - I(x). (7) implies [18, p.113] 

c(S%),£) = (£,x)+c{y l3 (0),£),\Jl£H (8) 

The definition of 7 _, (8) and c(^(0),£) = c(^(0),-£) 
imply 7_ = (£, x). Due to (4), x is the ^-estimation. 
Let us prove (5). If x £ ^(O) then 1^/3' ^x) < 1. Thus 
/3-5^(0) C Sf(0) implying ^(0) C /3^(0). If a: € 
S?(0) then Ii(/33a;) = fih(x) < 13 => /3^(0) C Sf^(0). 
Therefore 

8^(0) = /9*Sf(0) =► ctff p (p),£) = 0*c(Sf(O),O (9) 

Now (4) and (8) imply (5). Hence, £ = domc(Sf (0), •) 
due to Definition 1. 

Let us prove i?(T) =£.If^G #(T) then £ = L* z + H* u 
for some z £ £)(L*), u £ y and we get Vy £ ^(0) by 
Cauchy inequality [4, p. 4] 

(£, p) = (z, Lp) + (u, Hp) < \\z\\ 2 + \\u\\ 2 < +00 

so that R(T) C domc(^(0), ■) = £■ If £ $ R{T) then 
(£, x) > for some jc £ N(T*) as H = B(T) ® AT(T*). 
Noting that Sf(0) = {a; : ||T*a;|| 2 < 1} we derive 
c(&(0),£) > sup{(£,x) : T*x = 0} = +00. Hence, 
C c R{T). 

Let us prove (6). Set a := inf ve ^( a ) s&Vx&<s(-y) II V — x ll- 
Using Definition 1, one derives 

q = inf sup sup I (^, p> — x) I > 
sup inf sup I (£, y — x) I = sup p(£) = 

||i||=lV6*(»)xe«f(i/) ' ||£|| = 1 

sup sup |(£, x — x)| = sup ||x — x|| > a 

||i||=l x£W(y) xeW(y) 



Now assume £ = H. Since domc(^(0), •) = C, it follows 
c(Sf (0), •) is finite in % and therefore [9, §2.3] continuous. 
As a consequence, (6) is finite. 



4 7(2 -as) = 7i(s) + /(£)-2(Li,Lx) + 2(j/-iYx, J e r a;} 

5 Note [4, p.42] that \\p\\ = sup n=1 (£,(p) and [4, p. 55] 
inf x sup a F(x, y) > swp y mi x F(x, y) for convex-concave F. 



3 ^-estimation for non-causal DAEs 

Consider the model 

Fk+iXk+i - C k x k = fk+i,F x — /o, (10) 
Vk = H k x k + g k , k = 0, 1, . . . (11) 

where F k ,C k € M mxn , H k £ W xn , x k £ R n is a state, 
/ft G M" 1 is an input and y k ,g k £ IR P represent an output 
and the output's noise respectively. In what follows we 
assume that an initial state xq belongs to the affine set 
{x : F x = f }. We define £ T = ({/JoM&Jo) and 
assume 

T 

Hr £ <$ = Ur : ^{Sifi, fi) + {Ri9u9i) < 1} (12) 



PROOF. In order to apply Proposition 1, we rewrite (10)- 
(11) in the operator form: set U = (R") T+1 , y = 
(W>Y+\ T = (M m ) r+1 and p* = {x* a )l, y* = {y* s }l, 
f* = {f:Vo, V* = mi, 



F m „ 0„ 
-C F ± 0„ 

0™„ o"' o'' 



m „ m „ 

m „ m „ I ) JJ 

-C T _i F T 



Hq ... pn 



Define V T — (o»», ■ ••, o„„,b) and rewrite (12) with 
Qi = diag (So, ..., S T ), Q 2 = diag (R , ..., R T ) as 

ir = (/, tj) e Sf = {(/, r?) : (Qx/, /) + (Q 2 r?, r?) < 1}. 

It is clear that y* , H, g* , L, /*, p* , defined as above, 
satisfy (2)-(3) and (/*,»?*) € ^. Let #(j/*) denote the a 
posteriori set generated by y* . Then X(t) = V T (^(y*)) 
by definition. Thus 



where S* fc £ R mx ™ and i? fc € R pxp are positive definite 
self-adjoint matrices. 

Suppose we observe y{,. . . ,y* r , provided that y* k is de- 
rived from (11) with g k — g k and x k — x* k , which 
obeys (10) with f k = f*, and ({f*}h~, {g^) £ 0. De- 
note by X(t) the set of all possible states x T of (10) con- 
sistent with measurements y\ , . . . , y* and uncertainty 
description (12). 

Definition 2 We say that x T is a (.-estimation of the 
state x* in the direction £ £ M. n iff 

sup \(£,x — x T )\ = p(£,r) :— inf sup \(£, x — z)\ 

x£X(t) *G x (t)x€X(t) 

p(£, t) is said to be an £-error. A minimax observable 
subspace at the instant t for the model (10)-(11) is defined 
by C(t) — {£ : p(£, r) < oo}. I T :— n — dim. £(t) is called 
an index of non- causality of the model (10)-(11). 

Theorem 1 Define j3 T :— l — a T + (P T x T ,x T ) andx T = 
P+r T with r = H Q R y , a = (R yo,yo), 

Pk = H k R k H k + F k {S k _i — S k -iC k -iB k _ 1 C k _ 1 S k -i]F) 
Po = F SoF + H Q R Ho, B k = P k + C k S k C k , 
a k = a k -i + (R k y k ,y k ) - (B^rk-i^k-i), 
r k = FfcS , fc_iCfc_iBjJ'_ 1 rfc_i + H k R k y k , 

Then x T is the £-estimation of x*, C(t) — {£ : P^P T £ — 
£}, P (£,t) = J3?{P+e,Ql and 6 

X(t) = x t + J3?X(t),X(t) :={x: (P t x,x) < 1} (13) 



P (£,t)= inf sup \(£,V r (v-n=m (14) 

with / = V' r £. Hence, x T = V T (p, where (p is the 
Z-minimax estimation of the state p* of (3) in the 
sense of Definition 1. Proposition 1 implies <p = x 
and p(l) = /3h(&(0),l). Let us prove (13). (7) and 
(9) imply V T (&{v*)) = VrX + ^F T (&(Q)). Therefore, 
X(t) = V T x + /3^V T (^(0)). Now, let us prove x T = V T x 
by the direct calculation. Define 



V T (x , ...,X T ):= §(xo) + 22 $ «( x «' x 8 +i) 



with $ k (x,p) := \\F k+1 x - C fc p ||| fc + \\y k +i - 



H k+ix\\ Rk+1 , *(») :=||-Fba 



I S 



j/o - H x\\ 2 Ro , and 



6 Note that r >->• X(t) represents the a posteriori set-valued 
observer [20] 



B T (p) •■= min V t (xq . . .x T -!,p),B Q (p) := $(p) 

X ...X T _1 

Lemma 1. Let p £ W L . There exists (x\...x T -x) £ 
(R n ) T ~ 1 so that V T (x\ . . .x T -i,p) = B T (p) and 

B k (p) = (P kP , p) - 2 (r k , p) + a k ,P k >0 (15) 

Lemma 1 implies B T is a quadratic and non-negative 
function. Therefore AxgminS T = {x : P T x = r T } =/= 0. 
This and I(p) = V t (xq, ... ,x t ) imply 

B T (x T ) = minS T = min V T = min/ 

Defining x :— (x\ . . . x T —i,x T ) with x k taken as in 

Lemma 1 for p = x T , we obtain I(x) = mini and 

x T = V T x. Therefore, x T is a ^-estimation. 

We note mini = B T (x T ) = 1 - T . Thus (3 = J3 T by 

definition. 

Let us prove V T (&(0)) = X(t). Since ^(0) does not 

depend on y* and S?(0) = &(y*) provided y* = 0, 



we can calculate V T (^(0)) assuming y* — 0. In this 
case I — I\ — V r , x T = and B T (x) = (P T x,x) 
so that x € X(t) <& B T (x) < 1. If x € P T (&(0)) 
then, by definition, there exist x\ . . . x r ^\ so that 
V T {x\ . . .x T -i,x) < 1, implying B T (x) < 1. Now, let 
B T (x) < 1. Then V T (x\ . . .x T _i,x) — B T (x) due to 
Lemma 1 and thus x € V T (&{$)) by definition. 
Formulae (5), (14) and V T (&(0)) = X(t) imply 
p(£, T ) = 0c(X(t),£) = i3?{P+t,t)h for £ € C{t) and 
C{t) = domc(X(r),-) = {£ : P+P T i = £}. Details of 
calculation of c(X(t),-) are given in [18, p. 108]. This 
completes the proof. 

Proof of Lemma 1. We shall apply the dynamic pro- 
gramming [2]. Since V T is additive, it follows that 

B T (p) = min{$ T _i(x r -i,p) + S r _i(a; r _i)} (16) 

CCr-l 

V T is convex and non-negative by definition. Thus B T 
is non-negative and convex for any r € N. Convexity is 
implied by the definition of B T as for any convex func- 
tion (x,y) <-> f(x,y) the function y n- min^ f(x,y) is 
convex [18, p. 38]. 

Let us prove (15) by induction. (15) holds for B and P . 
We shall derive (15) for Bk+i, P k +i, assuming it holds 
for Bk, Pk- Define 

E k (x,p) := $ k (x,p) + (PkX,x) -2(r k ,x) + a k (17) 

Then B k+1 (p) = mm Xk E k (x k ,p) due to (16). Com- 
bining P k > with definition of <& k , we derive 
x h- > E k {x,p) is a convex quadratic function for any p. 
This and E k {x k ,p) > B k+1 {p) > imply [18, p.268] 
Argmin^Sfc^,^) ^ 0. On the other hand [18, T.27.4] 
x e Argmin a ,S fe (a;,p) iff V x E k (x,p) = 0. Finally, we 
obtain Argmin^ E, k (x,p) ^ and 

Argmin :E S fc (a;,p) = {x : B k x = G" k S k F k+1 p + r k } (18) 

If we set q k = B+{C' k S k F k+ ip + r k ) then q k £ 
Argmba. x Sk(x,p) due to [1]. Now, it is sufficient to 
calculate E k (q k ,p) in order to see that (15) holds for 
B k +i and P k +\. Assertion P k +i > holds since B k +i 
is convex. To conclude the proof, let us define x T = p 
and x k <E Argmm x E k (x,Xk+i), k — l,r— 1. Then 
V T (xi...x T -i,p) =B T (p) due to (16)-(18). 



3.1. Example. Consider a system p k +i = A k p k + 
v k ,po = v and assume y k = H k p k +g k provided p k G R n 
and (v, go, ... , g T ) belong to some ellipsoid. Now, given 
y* , . . . , y* , one needs to build the worst-case estima- 
tion of p* T . We cannot apply directly standard minimax 



framework [16,7,5] in this case as we do not have any 
information about the bounding set for (vq, . . . ,v T ). 
Instead, we apply the approach 8 , proposed above. De- 
fine F k := (E,0), C k ■- (A k ,E), H k := {H k ,0) and 
x* k := (p* k ,v* k ). Then x* k verifies (10)-(11) with f = v, 
f k = and g k — g k , k — l,r. Therefore, the original 
problem may be reformulated as: given y*, . . . , y*, one 
needs to build the I = (£, 0)-minimax estimation of x*. 
Of course, the estimation of x T in the direction I = (0, £) 
has an infinite minimax error for any £. But this is nat- 
ural as (I, x T ) = (v T ,£) and (vi, . . . , v T ) is unknown. 
In what follows we present a numerical example. 



Let p k G 



V k 



(" 



k sin(fc) 
10 



Hk — 

k sin(fc) 
10 



(i o), A k 



(I 



and 



). We have generated p* k with 



= (ii) and y* k with g k = 2 |^, R : - ; 



fe+i' 



S fc = diag(l,l), 

Fig. 1. 



k = 0, 50. The results are displayed on 




Fig. 1. The dashed line corresponds to the real values of 
(£,p k ) with I — (0, 1), k — 1,40; the solid line corresponds 
to the ^-estimation (£,Xk)\ the bold dashed lines represent 
dynamics of the boundary points of the segment £(X(r)). 
Note that the trajectory of the estimation is centered with 
respect to "the bounds" - bold dashed lines. 

3.2. Minimax estimator and Hijiix filters. In [5] 

a connection between set-membership state estimation 
and Hrx, approach is described for linear causal DAEs. 
The authors note that the notion of informational state 
(X(t) in our notation) is shown to be intrinsic for both 
approaches: mathematical relations between informa- 
tional states of -ffoo and set-membership state estima- 
tion are described in [5, Lemma 6.2.]. Comparisons of 
set-membership estimators with -/^/-Hoo filters for lin- 
ear DAEs are presented in [19], provided F k = E. Let 
us consider connections to ^-filters in details. In [10] 
the authors derive the Kalman's recursion to DAE from 
a deterministic least square fitting problem. Assuming 
rank [ H k ] = n, they prove that the optimal estimation 



7 x i-> (Ax, x) — 2(x,q) + c is convex iff A is a symmetric 
non-negative matrix. 



8 Since rank [ H k 
not applicable. 



< 2n, it follows that results of [10] are 



X%\k 



can be found from 
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Corollary 1 Let r$ = HqRqi/o- If ra,nk[ H k ] = n then 
4=0 and P^r k = x k \ k . 



PROOF. Let us set R k = E,S = E,S k = E for sim- 
plicity. The proof is by induction on k. For k = 0, P |o = 
Pq ■ The induction hypothesis is P k -i\k-i = P k \- Sup- 
pose A e W nxn , B e R nxn , B = B' > 0; then 



A(A'A + B' 



{E + ABA!)~ X AB (19) 



Using (19) we get ABA' = [E + ABA']A[A'A + 
B~ l ]~ l A 1 . Combining this with the induction assump- 
tion we get E + Cfc-iPk-iik-iCjU! = E + [E + 
Cfe-iPk_i|k_iC^,_ 1 ]Cfe_i[Pfe_i + C[_ 1 Ck_i] _ C' k 



fe-1- 
By simple calculation it follows from the previous 

equality that E - Ck_i(Pfc_i + qUiCfc-i) -1 ^-! = 
(E + Cfe-iPfc-iifc-iC^]^) -1 Using this and definitions 
of P k , P k \ k , we get P k x = P k \ k . 

Pq ro = xo|o due to corollary assumption. Sup- 
pose that P^^k-i = ifc— life— l ■ The induction hy- 



Pfcife and (19) imply (E 



pothesis, equality P k 
Cfe-iPfc-iifc-iC^j) - C k -ix k _\\ k _i = C k -i(C k _ x C k -i- 
Pk-i) k -i r k-i- Combining this with definitions of x k \ kl 
r k we obtain x k \ k = P,7 ffc- This concludes the proof. 



4 Conclusion 

We describe a set-membership state estimation ap- 
proach for a linear operator equation with uncertain 
disturbance restricted to belong to a convex bounded 
closed subset of abstract Hilbert space. It is based on 
the notion of an a posteriori set [16] &(y), informa- 
tional set [5] and the notion of the minimax observable 
subspace for the pair (L,H). The latter is new for the 
set-membership state estimation framework. It leads 
to nontrivial new results in set-membership state es- 
timation for linear non-causal DAEs: we present new 
equations describing the dynamics of the minimax re- 
cursive estimator for discrete-time non-causal DAEs. 
We prove that these equations arc consistent with the 
main results already established for regular DAEs. We 
illustrate benefits of considering non-causality in the 
state equation, applying our approach to a linear filtra- 
tion problem with unbounded noise. 
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